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A search for the Standard Model Higgs boson in proton-proton collisions at a center-of-mass energy of 7 TeV is 
presented in the decay channel H — > ZZ 2£2u. The search is conducted by the CMS experiment with data 
accumulated during the 2010 and part of the 2011 running periods of the LHC, for a total integrated luminosity 
of 1.1 fb~^. No excess is observed in the ZZ transverse mass. Limits are set on the production of the Higgs 
boson in the context of the Standard Model and in the presence of a sequential fourth family of fermions with 
high masses. 



^ 1. Introduction 

The search for the Higgs boson and its discovery are among the central goals of the experiments at the 
T-H LHC. Direct searches at the LEP collider have set a lower limit on the Standard Model (SM) Higgs mass of 

114 GeV/c^ [I], while direct searches at the Tevatron exclude the SM Higgs in the mass range 158-173 GeV/c^ 

^ at 95% confidence level [2]. The primary production mechanism for the Higgs at the LHC is gluon-gluon fusion, 

with a small contribution from vector boson fusion [T . 
JL In the following, a search for the SM Higgs boson in the H — > ZZ — > 2£2i/ channel is presented [3|. This 

channel is especially suited for high-mass Higgs boson searches, and a mass range from 250 to 600 GeV/c^ 
is considered in this analysis. Results are reported based on a data sample corresponding to an integrated 
'— ' luminosity of 1.1 fb~^, recorded by the CMS experiment at the LHC at a center-of-mass energy of 7 TeV. An 
^ interpretation of the results in terms of exclusion limits for the Higgs boson is also presented in a scenario 
^ with a fourth generation of fermions [S] added to the SM. For sufficiently large lepton and quark masses, this 
QQ SM extension has not been excluded by existing constraints (6j. Here we consider a fourth generation mass of 
Cn 600 GeV/c2. 
0^ 



2. The CMS Detector 



I The innermost component of the CMS detector is the tracking system, equipped with silicon pixel and 

• . microstrip detectors, used to measure the momenta of charged particles and reconstruct the interaction vertices. 

The tracker is immersed in a solenoidal magnetic field of 3.8 T and covers the pseudorapidity range \r]\ < 2.5. 
K> The tracker is surrounded by an electromagnetic calorimeter (ECAL), made of PbW04 crystals, and 

in ^ brass/scintillating fiber hadronic calorimeter (HCAL). Both calorimeters cover the pseudorapidity range 
C3 \t]\ < 3.0. To improve the detector hermeticity, forward quartz-fiber Cherenkov calorimeters (HE) extend 
the coverage up to \ri\ = 5.0. 

The outermost part of the CMS detector is the muon spectrometer, used to measure the momentum of muons 
traversing through the detector. It consists of three different types of gaseous detectors embedded in the iron 
return yoke: drift tube chambers (DT) in the barrel, cathode strip chambers (CSC) in the endcaps, and resistive 
plate chambers (RPC) in both regions. 

A detailed description of the CMS detector can be found in [7]. 



3. Simulation of Physics Processes 

In this analysis, signal and background datasets are produced by detailed Monte Carlo simulation of the 
detector response, based on GEANt4 [71 [5], taking into account the limited inter-calibration and alignment 
precision, and using the full CMS reconstruction chain. 

The H — >■ ZZ — 7> 2£2i/ signal and the Drell-Yan and tt backgrounds are generated using the program 
POWHEG [H]. The W + jets and single top events are simulated using madgraph [TO], while the diboson 
backgrounds {ZZ, WZ, WW) are simulated using pythia [TT]. In case of the irreducible ZZ background, a 
dynamic next-to-leading order fc-factor is applied as a function of the transverse momentum of the Z boson (see 
Eig. [l] left) . This fc-factor has been computed using the MCFM program [T^ . Contributions from the gg ZZ 
events are also taken into account [13] (see Eig.[l] right). 
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Figure 1: left: NLO to LO cross section ratio for ZZ — >■ 2£2i' at ^/s — 7 TeV as a function of the dilepton transverse 
momentum, computed using mcfm; right: contribution of the gg — ZZ subprocess relative to qq — )■ ZZ contribution. 



4. Event selection 

If the mass of the Higgs boson is much larger than twice the mass of the Z boson, each Z boson produced 
in the Higgs boson decay has a significant transverse momentum (pt)- Consequently, a iJ ^ ZZ — >■ 2^2i/ 
event is characterized by the presence of a boosted pair of isolated leptons, e+e^ or with invariant mass 

consistent with the Z mass, and large missing transverse energy (MET) arising from the decay of the other Z 
boson into neutrinos. 



4.1. Trigger and Preliminary Selections 

Since the final state involves two high-pT isolated leptons, events selected for the analysis are required to pass 
dilepton triggers. In case of the muon channel, a trigger with a threshold oi px > 7 GeV /c for each muon was 
used in the initial data-taking period (about 191 pb"""^), while for the remainder of the dataset the thresholds 
were 13 and 8 GeV/c for the leading and second muon, respectively. The efficiency for selecting signal events 
with these dimuon triggers is found to vary between 95 and 99%, depending on the transverse momentum and 
pseudorapidity of the muons. In case of the electron channel, the trigger has thresholds oi px > 17 GeV/c on 
the leading electron and px > 8 GeV /c on the second electron. The efficiency for this trigger is found to be 
larger than 99%. 

In order to suppress events with an incorrect MET measurement due to detector noise, a filter is applied to 
remove events in which anomalous HCAL noise is detected. Beam-halo events are also vetoed. 

An event is required to have at least one primary vertex within 24 cm of the geometrical center of the CMS 
detector in the beam direction and within 0.2 cm in the plane transverse to the beam. In events with multiple 
vertices, the one with the largest value of X]Pt fo'^ ^^e associated tracks is chosen as the reference vertex. 



4.2. Lepton and Dilepton Selection 

Muon candidates can be reconstructed using two algorithms: one uses tracks in the inner silicon tracker 
matched to segments in the muon system; the other algorithm uses tracks produced from a global fit of tracker 
and muon chamber hits, seeded by segments in the muon system fl4j . Further identification criteria based on 
the number of hits in the tracker and the muon system, the fit quality of the muon track and its consistency 
with the primary vertex, are imposed on the muon candidates to reduce fakes. 

Electron reconstruction also involves two algorithms, one in which energy clusters in the ECAL are matched 
to hits in the silicon tracker, and another in which tracks in the silicon tracker are matched to the ECAL 
clusters [TS]. These reconstructed electrons are required to pass further identification criteria, based on the 
ECAL shower shape, track-ECAL cluster matching, and consistency with the primary vertex. Additional 
requirements are imposed to remove electrons produced in photon conversions. 

To further suppress the QCD background, leptons are required to be isolated: a cone of radius Ai? = 
-f Arf = 0.3 is defined around the leptons, and the sum of transverse energy deposits in the calorimeters 
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and the transverse momenta of tracks in this cone is computed. This isolation smn is required to be smaller than 
15% (10%) of the momentum of the muon (electron). The median energy expected from pileup interactions is 
subtracted from the isolation sum, in order to reduce possible inefficiencies from enhanced detector activity in 
events with high pileup [16 . 

Events are selected such that there are two well-identified, isolated, opposite-charge leptons of the same flavor, 
each with pt > 20 GeV/c, which form an invariant mass consistent with the Z mass. In addition, events are 
required to have exactly two leptons with px > 10 GeV/c: events with any additional lepton are rejected in 
order to reduce the WZ background in which the W and Z bosons both decay leptonically. Moreover, the 
dilepton candidate is required to have a transverse momentum larger than 25 GeV/c. The impact of this cut on 
signal efficiency is small (less than 1% for all the Higgs boson masses considered), but allows the use of 7 + jets 
data to effectively model the Z + jets background, as explained in the following. 

With this selection, the principal backgrounds in this analysis are Z + jets, tt, single top, W + jets, WZ, 
WW, and ZZ events. 

4.3. B-Tag Veto 

In order to suppress backgrounds containing top quarks (single top, tt), events with at least one 6-tagged jet 
are vetoed. Particle Flow jets ^ with a transverse energy Et > 30 GeV are considered for 6-tagging. _B-jets 
are tagged using the "Track Counting High Efficiency" (TCHE) algorithm [TB], which uses displaced tracks in 
a jet to compute a 6-tagging discriminator. 

4.4. Missing Transverse Energy 

The signal process H — > ZZ — ^ 2121/ is characterized by a large missing Et, arising from the two neutrinos 
in the final state. Therefore, a high threshold is imposed on the Particle Flow MET 19j in order to suppress 
the large Drell-Yan background, which is characterized by little real missing transverse energy. However, jet 
mismeasurements and detector effects can produce large MET values in Z + jets events. Since this background 
is five orders of magnitude more abundant than the signal, the contamination of the selected signal region due 
to the high-MET tail is very significant. To further suppress backgrounds with fake MET coming from jet 
mismeasurements, events in which the MET is aligned with a jet are removed using a cut on the A(/)(MET,jet) 
variable. Both the MET threshold and the Aip cut are dependent on the Higgs mass hypothesis, since high-mass 
Higgs events are characterized by larger MET. 

4.5. Transverse Mass 

Higgs signal events show a narrower transverse mass distribution than background events, where the Higgs 
transverse mass is defined as 

Mt^ = (^^JpT.z^ + Mz^ + MET^ + Mz^^ - (pt,z + MEt) ^ . 

An m/f-dependent two-sided cut is thus applied to the Mt variable to further separate signal from background. 

The event selection criteria used in the analysis are hsted in Table |l] For MET, A(?!)(MET,jet) and Mr cuts, 
which are smoothly dependent on the Higgs mass, thresholds for selected mass points are shown as examples. 

5. Bacl^ground Estimation 

All the backgrounds in this analysis can be divided into three categories: 

• Z + jets with fake MET due to jet mismeasurement and detector effects; 

• non-resonant backgrounds (i.e., events without a Z resonance): tt, single top, WW, W + jets; 

• irreducible background: electroweak ZZ pair production and fully leptonic decays of WZ pairs. 

ZZ and WZ backgrounds are modeled using simulation, while the remaining backgrounds {Z + jets and all 
non-resonant ones) are estimated using data-driven methods. 
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Table I: Event selection cuts. 


Cut 


Cut Value 




PT > 20 GeV/c 


Z mass window 


|m« - 91.1876 GeV/c^l < 15 GeV/c" 


Transverse momentum of vetoed 3'^'* lepton 


PT > 10 GeV/c 


Zj LI Ojiio Vci be lliUlliCliLr U-lli 


Pt,z > 25 GeV/c 


6-tag veto (jet px > 30 GeV/c) 


TCHE discriminator < 2 


MET-jet separation 


mi/ = 250 GeV/c^ 
m/j = 350 GeV/c^ 
mff = 500 GeV/c^ 


A(^(MET,jet) > 0.62 rad 
A(^(MET,jet) > 0.14 rad 


Missing transverse energy (MET) 


ruH = 250 GeV/c^ 
= 350 GeV/c^ 
= 500 GeV/c^ 


MET > 69 GeV 
MET > 97 GeV 
MET > 141 GeV 


Transverse mass {Mt) 


= 250 GeV/c^ 
m^f = 350 GeV/c^ 
ruH = 500 GeV/c^ 


216 GeV/c^ > Mt > 272 GeV/c^ 
267 GeV/c^ > A/y > 386 GeV/c^ 
Mt > 336 GeV/c^ 



5.1. Z + jets Estimation 

The photon + jets process is very similar to Z + jets in terms of generating fake MET from mismeasured jets. 
Moreover, the cross section for 7 + jets events is much greater than that oi Z + jets. The former process can 
thus be used to build a model of the latter, provided that the two samples can be made kincmatically identical. 
This is achieved by rewcighting the 7 + jets events so that the pt spectrum of the photon matches the observed 
dilepton px spectrum. An additional reweighting is applied to match the jet multiplicity between the two 
samples. The photon is also assigned a random mass, according to the Z lineshape observed in dilepton events, 
and the yield of the photon sample is normalized to that of the Z. Lastly, the full analysis selection is applied 
to the reweighted photon sample. The MET distribution obtained with this procedure models very accurately 
the spectrum observed in Drell-Yan events (see Fig. [2|. The statistical error for the Z + jets prediction thus 
obtained is computed by summing in quadrature the weights of the 7 + jets events. 




Figure 2: MET spectrum of dimuon events (left) and dielectron events (right), compared witii reweighted single photon 
events in 1.1 fb~^ of data. 
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5.2. Non-Resonant Background Estimation 

The non-resonant backgrounds consist of several contributions. Events with and WW are characterized by 
a final state with two prompt leptons from W decays and real MET, plus two 6-jets for ti. Single top events 
produced in the most common tiy-channel show a similar signature. W + jets has instead one prompt lepton 
and a jet faking a second lepton. 

All non-resonant backgrounds can be estimated using events with a final state containing e+/i~ or e^/i+ pairs 
and passing the full analysis selection. The background level in the e'^e~ and final states is computed 

using the relations A^^^ = x N^.^ and N^f. = ae x Nf.^, where the scale factors a can be computed from the 
sidebands (SB) of the Z peak: 

]^SB njSB 
~ MSB ' ~ tsjSB ■ 

This method cannot distinguish between the non-resonant background and the signal in the H — ^ WW — > 212^ 
channel. Such signal events are therefore treated as a part of the background estimate. This signal contamination 
is found to be less than 7% of the total background yield for all considered Higgs masses. In the SM scenario 
with a fourth generation of fermions, this contamination increases, up to about 50% of the total background for 
a Higgs mass of 250 GeV/c^. Table O lists the predicted yields for non-resonant backgrounds with 1.1 fb~^. 



Table II: Yields from 1.1 fb"^ data for non-resonant 
Statistical uncertainties are also quoted. 





Predicted Yields 


MC Prediction 


250 


13 ±3 


12 ±0.65 


300 


3± 1.4 


4.8 ±0.39 


350 


2± 1.2 


1.4 ±0.21 


400 


0.68 ± 0.69 


0.42 ±0.11 


500 





0.038 ±0.033 


600 









:grounds in the muon (left) and electron (right) channels. 



mH 


Predicted Yields 


MC Prediction 


250 


9.4 ±2.2 


9.5 ±0.61 


300 


2.2 ± 1 


3.7 ±0.36 


350 


1.3 ±0.78 


0.99 ±0.16 


400 


0.5 ±0.5 


0.36 ±0.1 


500 





0.17 ±0.075 


600 









6. Systematic Uncertainties 

The signal and background yields are subject to several systematic uncertainties, summarized in Table |III| 
Most of these uncertainties apply to processes derived from simulation (Higgs signal, ZZ and WZ backgrounds) 
or to data-driven estimates of backgrounds (Z + jets and non-resonant backgrounds.) 



7. Results and Conclusions 

Figure [3] shows the Mt distributions of the selected dilepton events from data and the estimated background 
for the 250 GeV/c^ and 350 GeV/c^ Higgs mass hypotheses. No evidence of Higgs boson production is found. 
Table |IV| lists the event yields after the full selection for various Higgs mass hypotheses, for the estimated 
backgrounds, and for the 1.1 fb~^ dataset. 

Figure |4] (left) shows the 95% mean expected and observed C.L. upper limits on a x BR{H — > ZZ — )• 2l2v) 
[above), and the ratio R of such limit to the SM cross section csm (below), as functions of the Higgs mass mn- 
The results are obtained using a CL^ approach with a fiat prior for the cross section, for 1.1 fb^^. Such limits 
are not sufficient to exclude the existence of the SM Higgs boson in the examined mass range. The same results 
for the SM augmented with a fourth generation of fermions with masses of 600 GeV/c^ are shown in Fig. [4] 
(right). The presence of another fermion family produces an enhancement of the dominant gluon-gluon fusion 
cross section. With 1.1 fb^^, Higgs masses in the range 250 — 560 GeV/c^ can be excluded in this model. 

In this study, the exclusion limit is compared to the cross section for on-shell Higgs production and subsequent 
decay in the zero- width approximation. Acceptance estimates are obtained with Monte Carlo simulations based 
on ad-hoc Breit-Wigner distributions for describing the Higgs boson propagation. Recent analyses show that 
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Table III: Summary of the main systematic uncertainties affecting signal and background yields. 



Soiirc6 of iiric6rtciiiity 


XJncGrtciirity [/^] 


Luminosity 


6 


PDF, gluon-gluon initial state 


6-11 


PDF, quark-quark initial state 


3.3-7.6 


QCD scale, gluon-gluon initial state 


7.6-11 {ggH), 20 {ggZZ) 


QCD scale, quark-quark initial state 


0.2-2 (VBF), 5.8-8.5 {qqVV) 


Anti &-tagging 


1-1.2 


Lepton ID, isolation 


2 


Lepton momentum scale 


2 (2/x), 5 (2o) 


Jet energy scale 


1.5 


Pileup 


1-3 


Trigger 


2 (2^), 1 (2c) 


Non-resonant background estimation 


15-100 


Drell-Yan background estimation 


51-54 



Table IV: Expected number of signal and background events for an integrated luminosity of 1.1 fb~^ after full analysis 
selection, for dimuon {above) and dielectron {below) channels. Uncertainties for ZZ/WZ include the statistical and 
systematic components (in this order), non-resonant and Z-|-jets backgrounds have the full uncertainty quoted, and, 
finally, the signal column has only the systematic uncertainty quoted as the statistical uncertainty is negligible. 



mH 


zz 


wz 


Non- resonant 


Z + Jets 


Total 


Signal 


Data 


250 GeV/c^ 


6.3 ±0.1 ±0.5 


4.5 ±0.2 ±0.5 


13 ±3 


6.3 ±3.3 


30 ±5 


3.9 ±0.5 


28 


300 GeV/c^ 


5 ±0.1 ±0.4 


3.1 ±0.2 ±0.3 


3± 1.4 


3.5 ± 1.8 


15 ±2 


3.8 ±0.5 


13 


350 GeV/c^ 


5.2 ±0.1 ±0.5 


2.5 ±0.1 ±0.3 


2± 1.2 


3.5 ± 1.8 


13 ±2 


3.9 ±0.6 


11 


400 GeV/c^ 


3.6 ±0.1 ±0.3 


1.4 ±0.1 ±0.2 


0.68 ±0.68 


2.8 ± 1.5 


8.6 ± 1.8 


3.2 ±0.4 


6 


500 GeV/c^ 


2 ±0.9 ±0.2 


0.69 ±0.08 ±0.08 





1.6 ±0.9 


4.4 ± 1 


1.4 ±0.2 


3 


600 GeV/c^ 


1.1 ± 0.1 ± 0.1 


0.34 ±0.05 ±0.04 





1±0.5 


2.5 ±0.7 


0.57 ±0.08 


2 



niH 


ZZ 


WZ 


Non-resonant 


Z + Jets 


Total 


Signal 


Data 


250 GeV/c^ 


4.6 ±0.1 ±0.5 


3.6 ±0.2 ±0.4 


9.4 ± 2.2 


4± 2.1 


22 ±5 


3.1 ±0.4 


29 


300 GeV/c^ 


3.8 ±0.1 ±0.4 


2.2 ±0.1 ±0.3 


2.2 ± 1 


2.6 ± 1.3 


11 ± 3 


3.1 ±0.4 


21 


350 GeV/c^ 


4.2 ±0.1 ±0.4 


1.8 ±0.1 ±0.2 


1.3 ±0.8 


3.2 ± 1.7 


11 ± 3 


3.2 ±0.5 


8 


400 GeV/c^ 


3 ±0.1 ±0.3 


1.2 ±0.1 ±0.1 


0.5 ±0.5 


2.6 ± 1.3 


7.3 ±2 


2.5 ±0.3 


7 


500 GeV/c^ 


1.6 ±0.1 ±0.2 


0.62 ±0.07 ±0.07 





1.7±0.9 


3.9 ± 1.2 


1.1 ±0.2 


2 


600 GeV/c^ 


0.97 ±0.06 ±0.09 


0.31 ±0.05 ±0.04 





1±0.6 


2.3 ±0.7 


0.47 ±0.07 


1 



the use of a QFT-consistent Higgs propagator, and allowing for the ofF-shellness of the Higgs boson, dynamical 
QCD scales, and interference effects between Higgs signal and backgrounds, will result, at Higgs masses above 
300 GeV/c^, in a sizable effect on conventionally defined but theoretically consistent parameters (mass and 
width) that describe the propagation of an unstable Higgs boson [H |5D1 [5T]. In this work, these effects are 
included as an additional uncertainty of 10-30% (for 400 to 600 GeV/c^ Higgs mass) on the theoretical cross 
section. 
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Figure 3: Mt distribution of events passing the 250 GeV/c^ (above) and 350 GeV/c^ [below) Higgs mass selections, in 
the muon (left) and electron (nght) channels. 
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